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ABSTRACT 

The localization of the short-duration, hard-spectrum gamma-ray burst GRB 050509b by the 
Swift satellite was a watershed event. Never before had a member of this mysterious subclass of 
classic GRBs been rapidly and precisely positioned in a sky accessible to the bevy of ground-based 
follow-up facilities. Thanks to the nearly immediate relay of the GRB position by Swift, we began 
imaging the GRB field 8 minutes after the burst and continued for the following 8 days. Though 
the Swift X-ray Telescope (XRT) discovered an X-ray afterglow of GRB 050509b, the first ever of 
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a short-hard burst, no convincing optical/infrared candidate afterglow or supernova was found 
for the object. We present a re-analysis of the XRT afterglow and find an absolute position of 
R.A. = 12 h 36 m 13. s 59, Decl. = +28°59'04"9 (J2000), with a la uncertainty of 3'.' 68 in R.A., 3752 
in Decl.; this is about 4" to the west of the XRT position reported previously. Close to this 
position is a bright elliptical galaxy with redshift z = 0.2248 ± 0.0002, about 1' from the center 
of a rich cluster of galaxies. This cluster has detectable diffuse emission, with a temperature of 
kT = 5.25~t 16S keV. We also find several (~11) much fainter galaxies consistent with the XRT 
position from deep Keck imaging and have obtained Gemini spectra of several of these sources. 
Nevertheless we argue, based on positional coincidences, that the GRB and the bright elliptical 
are likely to be physically related. We thus have discovered evidence that supports the notion 
that at least some short-duration, hard-spectra GRBs are at cosmological distances. 

We also explore the connection of the properties of the burst and the afterglow, finding that 
GRB 050509b was underluminous in both of these relative to long-duration GRBs. However, we 
also demonstrate that the ratio of the blast-wave energy to the 7-ray energy is consistent with 
that of long-duration GRBs. This suggests a comparably high efficiency of 7-ray conversion as 
in long GRBs as might be expected if the same emission mechanism is at work in short and long 
GRBs. Based on this analysis, on the location of the GRB (40 ± 13 kpc from a bright galaxy), on 
the galaxy type (elliptical) , and the lack of a coincident supernova, we suggest that there is now 
observational consistency with the hypothesis that short-hard bursts arise during the merger of a 
compact binary (two neutron stars, or a neutron star and a black hole). In this context, we limit 
the properties of a Li-Paczynski "mini-supernova" that is predicted to arise on ~day timescales. 
Other progenitor models are still viable, and additional rapidly localized bursts from the Swift 
mission will undoubtedly help to further clarify the progenitor picture. 

Subject headings: gamma rays: bursts, gamma-ray bursts: individual: 050509b 

1. Introduction 

The distribution in duration (Mazets et al. 1981; Norris et al. 1984) and hardness (Kouveliotou et al. 
1993) reveals evidence for two distinct populations of classic gamma-ray bursts (GRBs): long-duration 
bursts, with typical durations around 30 s and peak energies at <~ 200 keV, and the minority short-duration 
bursts, with durations of a few hundred milliseconds (ms) and harder spectra. Despite remarkable progress 
in understanding the nature and progenitors of long-duration GRBs, comparatively little has been learned 
about the origin of short-hard bursts, primarily because very few such bursts have had rapid and precise 
localizations. 

The modeled bursting rate at redshift z = of long-soft bursts outnumbers short-hard bursts by about 
a factor of 3.5 in the BATSE catalog (Schmidt 2001); this assumes the same bursting rate as a function of 
redshift and does not include the effect of beaming, which, if different for long and short bursts, would imply 
that the intrinsic relative rates differ from those observed. While a number of bursts have been triangulated 
through the Interplanetary Network (see Hurley et al. 2005b) on roughly day-long timescales, there has only 
been one precisely localized short-hard burst relayed to ground observers in less than 1 hr (GRB 050202 / Swift: 
Tueller et al. 2005) 1 ; owing to its proximity to the Sun at time of localization, sparse groundbased followup 



There have been a few other short bursts (duration <2 sec) detected and well localized, but with soft spectra and hence 
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was undertaken. Including GRB 050509b, this corresponds to a ratio of 1:18 for short-hard to long-soft burst 
detections with Swift, much smaller than the BATSE result. 2 

As with long-duration bursts, the distribution of short bursts appears very nearly isotropic (Kouveliotou 
ct al. 1993; Briggs et al. 1996), and their brightness distribution (< V/V max >~ 0.35) is consistent with 
being a cosmological population. Still, there is no strong evidence to support the idea that short bursts are 
preferentially seen from z <; 0.37 rich Abell clusters (Hurley et al. 1997), nor are they clearly connected with 
star formation within ^100 Mpc (Nakar et al. 2005). 

Without precise and rapid localizations, the population statistics do not provide a strong constraint on 
the short-burst progenitors. Still, it has been largely reckoned that the leading candidates for short bursts 
are the merger of a neutron star binary (NS-NS; Blinnikov et al. 1984; Paczyhski 1986, 1991; Narayan et al. 
1992; Katz & Cancl 1996; Ruffert & Janka 1999; Rosswog & Ramirez-Ruiz 2002; Rosswog et al. 2003) or a 
black hole-neutron star binary (BH-NS; Lattimer & Schramm 1976; Eichler et al. 1989; Mochkovitch et al. 
1993; Kluzniak & Lee 1998; Bethe & Brown 1999; Popham ct al. 1999; Fryer et al. 1999). These systems hold 
several particular attractions. First, although uncertain, the estimated rate of mergers (between 1.5-20 per 
10 6 yr per galaxy; Belczynski et al. 2002; Sipior & Sigurdsson 2002; Rosswog et al. 2003) is comparable to the 
short-burst rate (Schmidt 2001). Second, the dynamical timescale of such mergers is several milliseconds and 
the sound-crossing times are of order ten milliseconds, comparable to the shortest observed bursts (Miller 
2005). Third, compact merger systems are likely to contain enough mass-energy in a transient torus to 
power short-burst flucnces as would be observed if at cosmological distances (Rosswog et al. 2003; Lee et al. 
2004; Rosswog 2005). The typical dynamical timescale in such binaries immediately prior to coalescence 
(ms) is much shorter than the observed burst duration, and so it requires the central engine to evolve into 
a configuration that is stable, while retaining a sufficient amount of energy to power the burst (Lee et al. 
2004). 

Mergers of such compact remnants are by no means the only possible channel to produce short bursts. 
Evaporating primordial black holes may produce short (< 100 ms) GRBs (Cline ct al. 1999), though basic 
energetics arguments suggest that it would be difficult to see such sources from distances well beyond the 
Galaxy. The recent discovery of a megaflare from SGR 1806—20 (Mcrcghctti ct al. 2005; Hurley et al. 2005a; 
Palmer et al. 2005; Terasawa et al. 2005) led to plausible suggestions that a substantial fraction (w40%) 
of short bursts could be produced by extragalactic magnetars (Hurley et al. 2005a). However, positional 
(Palmer et al. 2005; Nakar et al. 2005) and spectral (Lazzati ct al. 2005) arguments have led other workers to 
suggest that at most a few percent of the BATSE catalog could consist of short-burst magnetars. Note that 
not all compact mergers create fertile conditions (a transient torus around a BH) for making a short burst 
(e.g., Janka & Ruffert 2001; Rosswog et al. 2004). The duration of the burst in a compact binary merger is 
determined by the viscous timescale of the accreting gas, which is significantly longer than the dynamical 
timescale, thus accounting naturally for the large difference between the durations of bursts and their fast 
variability (Lee et al. 2004). In the collapsar scenario for long-duration bursts, on the other hand, the burst 
duration is given by the fall-back time of the gas (Woosley 1993; MacFadyen & Woosley 1999), which is 
typically greater than a few seconds. However, a modified collapsar scenario in which the burst duration is 
determined not by fall-back but rather by the dynamical timescales associated with the expanding outflow 



not members of the short-hard class. For example, GRB 040924 was a soft, X-ray rich GRB (Fcnimore et al. 2004; Huang et al. 
2005). Hereafter, we use the term "short burst" interchangeably with short-hard burst. 

2 As of 20 May 2005, Swift has localized 2 short bursts out of a total of 38; see 
http://swift.gsfc.nasa.gov/docs/swift/archive/grb_tablc.html . 
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might still meet the constraints of short GRBs (Woosley 2001). 

The theoretical predictions for the afterglows of short GRBs have been considered by Panaitescu et al. 
(2001). Since the peak flux of the prompt emission is comparable for short and long GRBs, if their distance 
scales are similar the isotropic equivalent energy output in gamma rays (-E 7i i so ) would be proportional to 
the duration of the GRB, which is ~ 10 — 100 times larger for long GRBs. If the efficiency for producing 
the gamma rays is comparable, then the isotropic equivalent kinetic energy in the afterglow shock (-Ek,iso) 
would have a similar ratio between long and short GRBs. This would imply the afterglow of short GRBs to 
be on average <~ 10 — 40 times dimmer in flux than that of long GRBs. The afterglows of short GRBs would 
be even much dimmer than this if they encounter a much smaller external density compared to long GRBs; 
this is the expectation from short-bursts from binary mergers outside of the host galaxy. Panaitescu et al. 
(2001) argued that a low external density would not affect the X-ray band, as the latter was assumed to 
lie above the cooling-break frequency, v c . We find that for a very low external density the electron cooling 
becomes very slow so that v c can lie above the X-ray band for the first few days, thus reducing the X-ray 
flux compared to that for a higher external density typical of the interstellar medium (ISM) found near 
star-forming regions of long-duration GRBs. 

To date the deepest early-time observations (At <; 1 hr) yielded upper limits Vjj m w 14 mag from the 
0.3 m ROTSE-I experiment (Kehoe et al. 2001). Hurley et al. (2002) compiled deeper non-detections at 
optical and radio wavelengths at times from days to weeks after four short bursts, with the faintest non- 
detection of R k, 22.3 mag at At — 20 hr (see also Gandolfi et al. 2000). Clearly, deep and early observations 
in search of a short-burst afterglow would require a rapid localization to an uncertainty comparable to the 
field of view of meter-class (and larger) telescopes. 

GRB 050509b (Gchrels et al. 2005) triggered the BAT coded-mask imager on-board Swift on 9 May 
2005 04:00:19.23 (UT dates and times are used throughout this paper; Hurkett et al. 2005). The position of 
GRB 050509b, with an uncertainty of 4' radius, was relayed to the ground within a few seconds. The initial 
localization was later revised to a position R.A. = 12 h 36 m 18 s , Deck = +28°59'28", with a 95% confidence 
error radius of 2.8' (Barthelmy et al. 2005a). Barthelmy et al. (2005a) describe the burst as a single-peaked 
source with duration of ~30 ms, peak flux of 2100 counts s _1 (15-350 keV), and a hardness ratio consistent 
with that of the short-hard population. At 06:29:23, a fading X-ray source was reported with a 6" localization 
(Kcnnea et al. 2005) and later updated to an 8" uncertainty radius at position R.A. = 12 h 36 m 13.9 s , Deck 
= +28°59'01" (Rol et al. 2005). 

GRB 050509b thus represents the first short-hard burst localized in real time to a position suitable 
for immediate follow-up observations from a suite of ground-based facilities. In this paper we describe the 
results of our observations of the field of GRB 050509b and what bearing these data have on the nature of 
short bursts and the physics of short-burst afterglows. In §2 we describe imaging and spectroscopy of the 
field. Our analysis of the X-ray afterglow of GRB 050509b is given in §3, leading to a localization near an 
elliptical galaxy (§4). In §5 we present a spectrum of that galaxy, its redshift, and inferred properties. We 
then argue, on statistical grounds, for a plausible association of this galaxy and the GRB. We demonstrate 
in §6 how GRB 050509b appears to be a subluminous burst relative to long-duration GRBs, but with a 
ratio of blast- wave energy to gamma-ray energy that is consistent with the long-duration population. In the 
remaining sections we describe new constraints on the nature of short-burst progenitors. Throughout, we 
assume a concordance cosmology with H — 70 km s _1 Mpc -1 , fl\ = 0.7, and Cl m = 0.3. All of the results 
presented herein, though generally consistent with our previous results in GCN Circulars, supersede them. 
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2. Observations and Reduction 

Initially, several groups reported (Rykoff et al. 2005; Ugarte et al. 2005; Bloom et al. 2005a; Torii 
2005) no new optical/infrared source that was consistent with the XRT position of GRB 050509b (Kennea 
et al. 2005). At 07:21:27 we highlighted the proximity of the XRT to a bright red galaxy (hereafter Gl = 
2MASX J12361286+2858580) and suggested a plausible physical association (Bloom et al. 2005a) based on 
its presumed membership in a z ~ 0.22 cluster (Barthelmy et al. 2005a). We later reported the determination 
of the redshift in Prochaska et al. (2005) and Prochaska et al. (2005). At 08:44:13 we noted the presence of 
a faint, compact source (hereafter SI; see 2) in the outskirts of Gl, which we deemed a plausible candidate 
counterpart (Bloom et al. 2005b). A very similar suggestion was made at 09:36:49 by Cenko et al. (2005c); 
in addition, they noted apparent variability of the candidate (later retracting the variability claim in Cenko 
et al. 2005b) and detection of three other faint sources (S2-S4) consistent with the XRT position (see also 
Cenko et al. 2005a) . Two additional sources (S5 and S6) in the XRT location were subsequently noted from 
Very Large Telescope (VLT) imaging by Hjorth et al. (2005), followed by another 5 sources (J1-J5) reported 
by Bloom et al. (2005). No radio emission (Parkinson 2005b; van der Horst et al. 2005) or GeV/TeV emission 
(Parkinson 2005a) is consistent with the XRT error localization. Below we discuss the observations, and 
further interpretation, leading to these reports. 

2.1. Optical and Infrared Imaging 

We observed the field of GRB 050509b on May 9 with the WIYN 3.5 m telescope and the OPTIC CCD 
imager with a 9.6'x 9.6' field of view and a plate scale of 0.14"/pixel. Under poor (~2") seeing conditions, 
two exposures totaling 360 s were obtained in the i' band beginning at 04.344 hr. In addition, we obtained 
2400 s of integration in the r' band under improved seeing conditions beginning at 06.088 hr. 

The data were reduced in the usual manner using flat-fields from both the illuminated dome and the 
twilight sky. The astrometric solutions to the individual images were calculated by comparison to the 
USNOB-1.0 catalog with a root-mean-square (rms) residual of 0.1". The photometric zero-points of the 
images were calculated by comparison to more than 50 stars in the Sloan Digital Sky Survey (SDSS) pho- 
tometry provided by Eisenstein et al. (2005). The zero-points of the WIYN images are uncertain at about the 
3% level. Limiting magnitudes were estimated from the histogram of fluxes in 10 4 seeing-matched apertures 
placed randomly within the field. The dispersion (a) of a Gaussian fitted to this distribution was used to 
estimate the 5(7 limiting flux in each image, which was converted to a magnitude using the known zero-point. 

The bright galaxy Gl to the west of the XRT position contaminates a significant portion of the 8" 
radius XRT error circle. We used galfit (Peng et al. 2002) to fit a smooth Sersic profile to this galaxy in 
order to remove most of the contaminant light prior to examining the XRT error circle. A series of 1000 
seeing-matched apertures placed randomly within the XRT error circle identified no new sources. The faint 
galaxy SI was detected at the > 5<r level in our deeper r' images. 

Near-infrared images were obtained with the 1.3 m PAIRITEL in the J,H, and K s bands (see Blake 
et al. 2005). Observations consisted of a 1130 s integration comprised of 7.8 s dithered exposures beginning 
at 04.1375 hr. These data were reduced by median-combining sets of individual exposures within a moving 
5-minute window. The resulting median was used to subtract the bright sky from the individual images. 
Finally, all of the individual images were combined to make high-resolution mosaics using a modified version 
of drizzle (Fruchter & Hook 1997). Zero-points were determined 2MASS stars in the field. Upper limits 
in the J,H,K S mosaics were estimated using the same technique as for the WIYN data. The WIYN and 
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PAIRITEL upper limits, as well as limits reported in the literature, are shown in Figure 1. 

We later imaged the field of GRB 050509b with the Keck I 10 m telescope and the LRIS-B instrument 
(Oke et al. 1995) using the dichroic D560 (50% transmission point at 5696 A) with G and R filters. Starting 
at 11.25 May 2005, beginning near astronomical twilight, we took 5 dithered images in each band for a total 
of 1660 s and 1620 s exposures in G and R, respectively. The data were reduced in the usual manner and 
combined, weighted by exposure times. 

On 17 May 2005, 08:05.5, 8.17 d after the GRB, we obtained deep i?(Ellis) (which is similar to Harris 
R; Bacon et al. 2003) imaging on the Echellete Spectrograph and Imager (ESI; Sheinis et al. 2002) on the 
Keck II 10 m telescope. In the presence of bright glare from the Moon, we combined several reduced images 
for an effective exposure time of 960 s. Since there is a negligibly small color term in converting Harris R to 
i? c , 3 we found a zero-point relative to the LRIS R image. There are no new sources to R c « 25.0 mag (5a), 
nor significant variations of the faint sources in the XRT error circle. 

3. The X-ray Emission 

The Swift XRT (Burrows et al. 2000) began observations of GRB 050509b on 2005 May 9 at 04:00:56, 
approximately 61 s after the BAT trigger. The observations consisted of eleven blocks, each about 2.5 ks in 
duration (except the first observation of 1.6 ks and the last observation of 1.8 ks), spread over a period of 
~21 hr. The XRT operated in a number of different modes throughout the observations. The most common 
(32.3 ks of exposure) and most useful mode for this object was the "Photon Counting" mode, which retains 
the full imaging and spectroscopic resolution of the instrument. The images are 480 x 480 pixels, with a 
scale of 2'.'36 per pixel. The XRT point-spread function is energy dependent, with a half-power diameter of 
18" at 1.5 keV. The energy resolution is also a function of energy, varying from about 50 eV at 0.1 keV to 
about 190 eV at 10 keV. 

The first Photon Counting observation began at 04:01:20 (Kennea et al. 2005) and lasted 1640 s. As 
noted in Kennea et al. (2005) and Rol et al. (2005), a faint X-ray source is detected in this first of the eleven 
observations, but it faded quickly below the background. We have obtained the XRT data from the Swift 
archive, and have analyzed them to determine the position of this X-ray afterglow candidate as well as to 
examine its variability. We briefly review the data reduction, and then we discuss the localization of the 
afterglow candidate and attempt to quantify the decay. 

3.1. Swift Data Reduction 

Using the Level 1 data from the Swift archive, we ran the xrtpipeline script packaged with the HEAsoft 
6.0 software supplied by the NASA High Energy Astrophysics Science Archive Research Center 4 . We used 
the default grade selection (grades to 12) and screening parameters to produce a Level 2 event file re- 
calibrated according to the most current (as of 2005 May 15) calibration files in the Swift database 5 . To 
produce images for source detection, we used the xselect software (also part of HEAsoft 6.0), with a filter to 

3 http://www. ast.cam.ac.uk/~wfcsur/technical/photom/colours/ . 
4 http://heasarc. gsfc.nasa.gov/ . 

5 http://heasarc.gsfc. nasa.gov/docs/heasarc/caldb/swift/ . 
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include only counts in PI channels 30-1000 (corresponding to photon energies of 0.3-10 keV). The PI channel 
to photon energy conversion was accomplished with the redistribution file swxpc0tol2_20010101v007.rmf 
from the calibration database. The effective area of the XRT at the position of the afterglow candidate was 
determined with the xrtmkarf tool, using the correction for a point source. 

3.2. X-ray Afterglow Localization 

A number of factors make the localization of this X-ray afterglow difficult. It is intrinsically faint and 
superposed on diffuse X-ray emission from a galaxy cluster at z = 0.22 (Gal et al. 2003). The initial source 
detection was performed with the wavelet-based routine wavdetect (Freeman et al. 2002), supplied with the 
CIAO 3.2 software package, which in our experience is quite good at detecting faint sources. We chose 
parameters appropriate for detecting point sources in this XRT observation; the pixel scales considered were 
a \[2 series starting at 4 pixels (4, 5.657, 8, 11.314, 16), and the significance threshold was set at 4 x 10~ 6 , 
corresponding to a ^1 false positive detection of a point source in the image. We detect 22 compact sources 
in the entire 32.3 ks data set. 

To study the properties of the afterglow candidate, we extracted the all the events within an area of 
radius 10 pixels around the nominal wavdetect position. In the first observation of 1.6 ks, there are 14 
counts in this region. When examining a plot of the cumulative distribution versus time, we noticed that 
the majority of the counts from this region occurred in the first 300 s. We therefore further investigated this 
brief interval. 

In the first 300 s of the first Photon Counting observation, the XRT detected 92 counts on the entire 
chip, with 73 of them outside of the 22 source regions. Within any 10-pixcl radius source region, we therefore 
expect an average of 0.1 background counts. We detect 9 counts in this region of the X-ray afterglow, with 
a reasonable expectation that all 9 are from the X-ray afterglow. Using the mean location of just these 9 
counts, we can obtain a relatively uncontaminated estimate of the source position. We calculate the 68% 
confidence interval in each direction as Tcrj/VA, where N is the number of counts (9), <jj is the sample 
standard deviation of the 9 coordinates in each direction, and [— T,T] is the 68% confidence interval of the 
Student's t distribution with N — 1 degrees of freedom. This gives us a position estimate, in the Swift XRT 
reference frame, of R.A. = 12 h 36 m 13 s 94, Deck = +28°59'05'.'3 (J2000) with an uncertainty of 3'.'6 in R.A. 
and 3'.'5 in Deck This is 4'.'3 North of the revised XRT position reported by Rol et al. (2005). A possible 
reason for this offset is that the Rol et al. position is based on 6.6 ks of XRT exposure and thus includes 
contributions from the diffuse cluster emission (see Figure 5), biasing the position estimate. 

We examine the absolute astrometric accuracy of the Swift XRT frame by searching for possible coun- 
terparts of the other 21 XRT sources in deep optical images. The best suited optical data for this is a Bok 
S-band image (Engelbracht & Eiscnstein 2005) because it covers an area large enough to contain the entire 
XRT field. Using a cross-correlation to 250 2MASS positions, we fit an absolute WCS using IRAF/CCMAP 6 . 
The overall geometry plus the considerable distortion across the Bok £?-band image was well fit by a fourth- 
order polynomial with rms residuals of 0.135" in R.A. and 0.158" in Deck Assuming a 100 mas global 
uncertainty in the 2MASS-Intcrnational Celestial Reference System (ICRS) tie, the absolute astrometry in 
the wide-field optical frame is thus uncertain to 170 mas in R.A. and 187 mas in Deck 



6 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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For each XRT source other than the afterglow candidate, Figure 3 plots the offset between the XRT 
position and the position of the closest optical source. Two XRT sources had two optical sources within 5"; 
for these, the closest optical source is represented by dashed lines and the next closest by dotted lines. There 
is an obvious locus around a 4". 5 difference in R.A., suggesting that these XRT sources are associated with 
the corresponding nearest optical sources. At a detection sensitivity around 10 -14 erg cm~ 2 s _1 , it is not 
surprising to find so many optical counterparts in the moderately deep Bok image. In a Chandra/ 'Subaru 
study of the R.A. = 13 hr XMMj R OS A T field, McHardy et al. (2003) find unambiguous optical counterparts 
for 61 of the 66 X-ray sources above 10~ 14 erg cm~ 2 s _1 . The mean R magnitude of these sources is R — 20.7, 
and the faintest counterpart is at R = 24.4 mag. 

Using the 14 sources in the above locus (excluding the two sources with multiple possible counterparts), 
we derive an offset between the XRT frame to the optical frame of 4'.'49 ± 0'.'72 W in R.A. and 0'.'42 ± 
0'.'30 S in Decl. Our best estimate for the location of the X-ray afterglow is therefore R.A. = 12 h 36 m 13 s 59, 
Decl. = +28°59'04"9 (J2000); this is 4.1" west and 3.9" north of the revised XRT position reported in Rol 
et al. (2005). The uncertainty in our position is a combination of the statistical uncertainty of the XRT 
localization (3"6 in R.A., 3'.'5 in Decl.) and the uncertainty in shifting the XRT frame to the ICRS (0"76 in 
R.A., 0'.'40 in Decl.). 

The astrometry in our original reports from WIYN and Keck imaging were based on a frame of approx- 
imately 10 stars in the 2MASS catalog. The release of the SDSS data and calibrations of this field allow us 
to improve the astromctric tie to the ICRS. We fit the Keck/LRIS G-band image to 91 sources in common 
with the SDSS object catalog with a third-order polynomial solution using IRAF/CCMAP. The uncertainty 
in the astrometric tie to SDSS, based upon residuals from the fit, is <r(R.A.) = 0.134"and cr(Decl.) = 0.153". 
Assuming a 75 mas astrometric uncertainty in the SDSS astrometric calibration to the ICRS (Pier et al. 
2003), we estimate the absolute uncertainty in the Keck-ICRS tie is <r(R.A.) = 0.154"and (j(Decl.) = 0.171". 

The XRT location is 11.2"± 3.6" (or 40 ± 13 kpc in projection) from Gl as we first noted in Bloom 
et al. (2005a). Spectroscopy of this source reveals that it is indeed an early-type galaxy (see § 5) and is a 
member of a cluster NSC J123610+285901 at z w 0.22 (Gal et al. 2003; Barthclmy et al. 2005a). Near the 
location of the revised XRT error circle, we find ^11 faint sources (all of which we or others have reported 
previously; see above). Figure 2 shows the Keck G and R images with identified source labeled. Table 2 
gives the astrometric positions and magnitudes of the sources. 

3.3. X-ray Afterglow Decay 

We examine the first 1.6 ks block of observations to characterize the temporal properties of the X-ray 
afterglow. A Kolmogorov-Smirnoff (K-S) test on the arrival times of the 14 photons gives a probability of 
0.06% that they come from a source with constant count rate. The next step in model complexity is one in 
which the X-ray count rate Rx in this region has a constant component (due to the background and diffuse 
cluster emission) plus a component with a power-law dependence on time (due to the fading afterglow). Our 
model is thus Rx{t) — A(t — to)~ a + B, where B is the constant (background plus cluster) count rate, to is 
the time of the BAT trigger, and A is a normalization chosen such that the model preserves the detected flux 
over the 1.6 ks under consideration. We determine B to be 0.00107 count s _1 from the later observations. 

We considered a range of a from to 4 and computed the K-S probability of the observed data coming 
from the model for each value of a. The K-S probability was highest (97.8%) at a w 1.3. For a < 1 
(0.77) and a > 1.7 (2.1), the K-S probability dropped below 32% (5%). For a = 1.3, the normalization 
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A is 22 count s" 1 . We can translate this to an energy- flux normalization by determining the conversion 
from counts to erg cm" 2 . We consider only the first 300 s of data for this determination in order to reduce 
contamination from the background. For each of the 9 counts, we know its energy as well as the effective 
area of the XRT at that energy. The average is 3.14 x 10" 11 erg cm" 2 count -1 . Our model for the X-ray 
flux (0.3-10 keV) of the afterglow only is then Fx(t) = A'{t — to)~ a . For a = 1.3, the normalization is 
A' = 6.9 x 10" 10 erg cm" 2 s" 1 . For example, the X-ray flux of the afterglow at t = 200 s after the BAT 
trigger for a = 1, 1.3, and 1.7 is Fx{200) = 5.8, 7.0, and 6.5 x 10" 13 erg cm" 2 s" 1 , respectively. 

Figure 4 shows, on a common scale, X-ray light curves for a number of GRBs and core-collapse super- 
novae. The location of the X-ray transient associated with GRB 050509b, when placed at z = 0.2248, is 
striking. While the slope of the transient agrees well with those of typical GRBs, its flux falls well below 
the typical long-GRB range. In fact, for any reasonable redshift (i.e., z < 3-5), GRB 050509b would still 
be significantly underluminous in its X-ray afterglow when compared to those of long-duration GRBs (see 
Figure 4). For the assumed redshift of the tentative host galaxy, the extrapolated X-ray luminosity at a 
few days, which is also consistent with the Chandra upper limit (Patel et al. 2005), is close to those seen in 
typical core-collapse supernovae. 

3.4. Diffuse Galaxy Cluster Emission 

We used wavdetect to search for large-scale structures in the full 32.3 ks XRT data set. The pixel scales 
searched were (20, 28.28, 40, 56.57, 80). The center of the diffuse emission presumably associated with 
the galaxy cluster had a wavdetect-determined position of 12 h 36 m 18 s 26, +28°59'06'/7. Figure 5 shows an 
adaptivcly smoothed image (using the CIAO tool csmooth) of the XRT data with the cluster center and 
GRB indicated. The colors represent the 0.3-10 keV count density. Contours are drawn at 0.00449, 0.00646, 
0.00934, 0.0136, 0.0197, and 0.0273 count arcsec 2 . As the image shows, the wavdetect-datGiuiniGd position 
of the diffuse emission is about 14" to the west and 4" south of the peak of the diffuse emission, which is at 
12 h 36 m 19 s 33, +28°59'10'/8 (J2000). [Note that the optical cluster center is 12 h 36 m 10 +28°59'00.9" (J2000) 
as defined by the center of the galaxy overdensity; this is about 125" east and 10" south of the peak of the 
diffuse X-ray emission.] We thus find that the XRT afterglow position is 75" west, 6" north of the cluster 
center, as defined by the peak of the diffuse X-ray emission, about 270 kpc in projection. 

We extract a spectrum from a region of 110" in radius centered on the wavdetect position. We use 
a similar-sized region in a source-free area to extract a spectrum for background subtraction. We require 
the cluster spectrum to contain at least 20 counts per bin, and we consider the range 0.3-10 keV. We fit 
the background-subtracted cluster spectrum in Xspec vl2.2 (Arnaud 1996) with a MEKAL (warm plasma) 
model absorbed by a Galactic column density of 1.52 x 10 20 cm -2 (Dickey & Lockman 1990). We set the 
MEKAL redshift at z = 0.2248 and the metallicity at [Fe/H] = 0.26 (Mushotzky & Loewenstein 1997) and 
allow the temperature and normalization to vary. The best fit temperature is kT = 5.25t 3 ;gg keV, which 
gives x 2 /dof = 22.4/20. 

4. Associating GRB 050509b with Gl 

We are now in a position to explore the possible association of GRB 050509b with the cluster and with 
the nearby elliptical galaxy Gl. Focusing on the BAT localization alone, we first consider the probability 
that a random position in the sky would be in a rich cluster of galaxies (here we neglect the effects of lensing, 
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expected to be small; for example, Grossman & Nowak 1994). A reasonable estimate of the covering fraction 
of clusters on the sky is given by the the DPOSS Northern Sky Optical Cluster Survey (Gal et al. 2003). 
Although this survey is not very deep (zn m <~ 0.3), low-redshift clusters should dominate the sky density. 
Gal et al. (2003) find a covering fraction of ~ 0.03 assuming a typical cluster radius of 1 Mpc, which suggests 
a chance alignment is improbable but not impossible. Moreover, the XRT localization of GRB 050509b to 
within 45" of the center of such a cluster would occur by chance with a probability of just ~7x 10~ 4 , but 
it is difficult to estimate a posteriori how large a distance from a cluster center one would have considered 
"significant." 

While the gas from the cluster environment may enhance the probability of localizing short-burst after- 
glows (see below), our expectation is that short GRB progenitors are caused by the death of stars of some 
sort, with the burst rate determined by processes on scales significantly smaller than cluster lengths. To this 
end, we should consider the chance probability of the GRB event occurring at close impact parameter to a 
galaxy similar to Gl. As reported by Eisenstein et al. (2005), the galaxy Gl has a Petrosian r' magnitude 
of 17.18 ± 0.02 mag based on imaging by the SDSS. The sky density of galaxies with comparable apparent 
magnitude brighter than Gl is <~ 40 per square degree (Blanton et al. 2003b). Therefore, the probability of 
an event randomly occurring within 20" (about twice the observed offset) of this bright galaxy is ~ 5 x 10~ 3 . 
We consider this a conservative estimate because this probability makes no reference to the galaxy redshift, 
type, size, or age (which are consistent with a priori discussions of short-burst host galaxies; e.g., Bagot 
et al. 1998 and Bloom et al. 1999). 

If one argues that GRB 050509b is indeed physically associated with this bright, low-redshift elliptical 
galaxy, one must consider why the several other well-localized short bursts have not shown similar associ- 
ations. The first possibility, that the short bursts arise from a more local population (as suggested by the 
magnetar flare from 27 December 2004; Hurley et al. 2005a) and GRB 050509b must therefore arise from a 
different population, was discounted for four of the best-localized short bursts (Nakar et al. 2005). Another 
possibility is that GRB 050509b was significantly closer than the other well-localized short bursts. A strong 
test of this hypothesis is to determine if other short bursts are associated with more distant clusters or in- 
trinsically bright, massive galaxies (e.g., through a deep imaging campaign). The third possibility is simply 
that short-burst progenitors need not always arise in such galaxies. In fact, for the NS-NS hypothesis we 
would expect mergers in galaxies spanning a wide range of Hubble types. A delayed BH-NS merger is also 
possible, but less likely if GRB 050509b is associated with Gl: statistically, the distribution of timescales 
for BH-NS coalescence is as broad as that for NS-NS coalescence, albeit quite model dependent (Belczynski 
et al. 2002; Sipior & Sigurdsson 2002), but the systemic kick velocity is expected to be systematically lower 
by a factor of a few in most theoretical models of formation of BH-NS binaries (kick velocity is roughly 
inversely proportional to mass, so more massive binaries receive less kick). Moreover, larger velocity kicks 
generally lead to shorter merger times. Thus the ~40 kpc offset tends to favor NS-NS over BH NS mergers. 

While these possible associations are tantalizing, a posteriori statistics are very suspect. Had the 
GRB been near a bright spiral galaxy, we might have made similar claims based on chance probabilities. 
Nevertheless, it remains the case that many workers had predicted the distinct possibility that a well-localized 
compact merger and/or short burst could be near an elliptical galaxy (Bagot et al. 1998; Bloom et al. 1999; 
Panchenko et al. 1999) (see also Livio et al. 1998), and so we suggest that these arguments might reasonably 
reflect a true association. Moreover, the possible association with an early-type host stands in stark contrast 
to results from long-duration GRBs (Bloom et al. 2002; Le Floc'h et al. 2003; Djorgovski et al. 2003), and is 
reminiscent (e.g., Bloom et al. 2002; Dar 2004) of the dichotomy between core-collapse and thermonuclear 
(Type la) supernovae. Only a larger sample of short GRBs will provide truly compelling evidence for such 
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a parallel. Still, based on the arguments above, we proceed by accepting the hypothesis that GRB 050509b 
is physically associated with Gl. 



5. A Putative Host Galaxy at z=0.2248 

We obtained a spectrum of Gl with DEIMOS (Faber et al. 2003) on the Keck II 10 m telescope under 
photometric conditions. The data were acquired in a series of two exposures starting at 07.47 hr on the night 
of the burst. The instrumental setup included the 600 line mm -1 grating blazed at 7500 A and centered 
at 7200 A, the GG455 order-blocking filter, and standard CCD binning. This setup gives nearly continuous 
wavelength coverage in the range 4500-9000 A. We observed the galaxy though a 1.1" x 20" slit at sky 
position angle 90° and an airmass of 1.0. This setup yields a FWHM resolution of <~ 5 A (i.e., a = 100 
km s _1 ). The data were reduced and calibrated with the DEEP spectroscopic pipeline for DEIMOS data 
(Cooper et al. 2006). Wavelength calibration and flat-fielding were performed using spectra of Xe-Ne-Kr-Ar 
and quartz lamps (respectively) obtained that night. 

The software provides a two-dimensional, sky-subtracted image of the spectrum across two CCDs of the 
DEIMOS mosaic. Unfortunately, the CCD that includes the bluest data has a pair of blocked columns which 
lie near the center of the galaxy profile. Therefore, we extracted the one-dimensional (ID) spectrum on this 
CCD using optimal extraction techniques assuming a Gaussian profile with a = 9.2 pixels (i.e., 1.1"). For the 
other CCD, we extracted a ID spectrum by adopting a 26-pixel (3") boxcar aperture. Finally, we processed 
and calibrated a spectrophotometric standard star (BD+28°4211) observed at the end of this night. After 
comparing its observed flux (in digital numbers) against the STIS CALSPEC calibration 7 , we calculated a 
sensitivity function which could be applied to our galaxy spectra. 

Spectroscopic observations of Gl, SI, S2, and 2 unidentified sources were obtained with the GMOS 
spectrometer (Hook et al. 2004) on the Gemini North 8 m telescope beginning at 10.27 May 2005 under 
photometric conditions. We used a 0.75" slit, a R400 grating blazed at 7640 A, GG455 order-blocking filter, 
and set the central wavelength to 6500 A. The airmass was low (1.0-1.1), so the effects of atmospheric 
dispersion were negligible (Filippenko 1982). Standard CCD processing and spectrum extraction were ac- 
complished with IRAF using a 1.74" aperture for these sources of interest (S1,S2). The data were extracted 
using the optimal algorithm of Home (1986). Low-order polynomial fits to calibration-lamp spectra were 
used to establish the wavelength scale. Small adjustments derived from night-sky lines in the object frames 
were applied. Using techniques discussed in Wade & Home (1988) and Matheson et al. (2000), we employed 
IRAF and our own IDL routines to flux-calibrate the data and to remove telluric lines using the well-exposed 
continua of the spectrophotometric standard EG-131 (Bessell 1999). 

Figure 6 presents the ID flux-calibrated spectrum of Gl against a vacuum, heliocentric-corrected wave- 
length array. The dotted line traces a la error array based on Poisson counting statistics. We have marked 
a number of detected absorption-line features and also the expected position for several strong transitions 
frequently observed in emission-line galaxies (e.g., Ha, [O III]). We have fit a double-Gaussian profile to 
Ca II H&K and measure z — 0.2248 ± 0.0002. This is consistent with the redshift inferred photometrically 
for this cluster from DPOSS (Gal et al. 2003). At this redshift, the luminosity distance is 1117.4 Mpc, and 
1" corresponds to 3.61 kpc in projection. 



f tp://f tp. stsci .edu/cdbs/cdbs2/calspec/ 
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To estimate the velocity dispersion of the galaxy, we have compared the spectrum against a template 
spectrum of HD 72324 (e.g., Kelson et al. 2000) smoothed by a wide range of a. The best match to the 
absorption lines of Gl with A reS i = 4000-5300 A is a = 275 ± 40 km s _1 . Accounting for the instrumental 
resolution, we derive a light-weighted velocity dispersion for this galaxy of 260 ± 40km s _1 . 

The spectral features evident in Figure 6 are typical of early-type galaxies. The spectral type and velocity 
dispersion indicate a massive elliptical galaxy with no apparent ongoing star formation. A quantitative limit 
to the current star-formation rate (SFR) can be inferred from the upper limit to the Ha luminosity of this 
galaxy. The emission-line flux in a 10 A window (Av ~ 300 km s _1 ) centered at the expected wavelength 
of Ha has a 3a upper limit of 1.2 x 10~ 16 erg s _1 cm~ 2 . Adopting the current concordance cosmology, we 
derive an Ha luminosity Lu a < 1.2 x 10 40 erg s _1 . Using the empirical relation between SFR and Lua 
(Kennicut 1998), the 3tr upper limit to the current SFR is 0.1 Mq yr _1 . 

A morphological fit to WIYN /-band imaging using galfit (Peng et al. 2002) shows good agreement with 
a de Vaucoulcurs profile (Sersic index = 4), with a x 2 /dof = 1.22. The effective radius is R e = 0.96" = 
3.47 kpc. The galaxy has an axis ratio of 0.81 with the semimajor axis aligned along a position angle east of 
north at ~90°. There was little improvement in x 2 /dof by adding more complicated morphologies or letting 
the Sersic index vary. 

The coincidence of a point source at radio wavelengths with the optical center of Gl might suggest the 
presence of a low-level active galactic nucleus despite the lack of telltale features observed in the Gl optical 
spectra. Moreover, inspection of archival images of this galaxy from the Near-Earth Asteroid Tracking 
Program (Pravdo et al. 1999) on 9 April 2002, 20 April 2002, 3 May 2002, 22 March 2003, and 8 April 2003 
reveals no apparent variability of the optical light from Gl. However, radio emission without corresponding 
optical emission is not uncommon in giant elliptical galaxies harboring mildly active nuclei (Ho 1999). The 
radio emission in Gl is unlikely to be associated with star formation, given the low SFR deduced above. 

The properties of this probable host galaxy contrast significantly with those measured for the galaxy 
hosts of long-duration GRBs. First, most hosts of long-duration GRBs exhibit emission- line features in- 
dicative of high SFRs (e.g., Djorgovski et al. 2003). Second, the absolute if -band luminosity of this galaxy 
(«1.6 xlO 11 /©) exceeds that of all previously identified GRB host galaxies (Chary et al. 2002). Third, the 
impact parameter of the GRB (as defined by the 90% XRT error circle) is larger than that of all previously 
associated GRB-host galaxy pairs (long-burst offsets <;10 kpc; Bloom et al. 2002). 

5.1. SI, S2: Faint Blue Galaxies in a High-Redshift Group? 

The Gemini/GMOS spectra of SI and S2 are featureless and blue. Examining the regions of the spectrum 
where Ha or R(3 would lie if at the redshift of the cluster, we detect no measurable emission. Assuming, 
for the moment, that the sources are at the cluster redshift of 0.22, we put a 3cr upper limit on the Ha 
luminosity of Ln a < 1.5 x 10 39 erg s _1 and Ln a < 1.4 x 10 39 erg s _1 for SI and S2, respectively. Using 
the equation from Kennicut (1998) relating the Ha luminosity to the SFR, we find that the upper limits for 
the unextinguished SFR, assuming that SI and S2 are cluster members, are ~ 1.1 x lO~ 2 M yr -1 for the 
galaxies. If SI and S2 are cluster members, then they are not forming stars, which would seem to conflict 
with their blue colors. 

A more likely scenario, also mentioned by Cenko et al. (2005a), is that SI and S2 are both background 
galaxies. Although the Gemini spectra range from 4600 A to 8600 A, the data have poor signal-to-noise ratio 
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blueward of 5200 A. Nevertheless, the spectral slope is well constrained and it suggests that these galaxies 
are forming stars (i.e., the continuum is relatively blue). The lack of corresponding emission lines (Ha, H/3, 
[O III] A5007, [O II] AA3727) falling in our spectral window therefore suggests that SI and S2 have z J> 1.3. 
Additional spectroscopy will be required to confirm our hypothesis that SI and S2 are faint blue galaxy 
members of a small group at moderate redshift. 



6. Theoretical Interpretation 

The fluence of the prompt gamma-ray emission measured by the Swift BAT is / = (2.3 ± 0.9) x 
10~ 8 erg cm~ 2 (Barthelmy et al. 2005b), which at the redshift of the tentative host implies an isotropic 
equivalent energy output of -E 7 .iso = (2.7 ± 1) x 10 48 erg. Since vF v is still rising roughly as v 5 in the 
15-150 keV Swift range, the total fluence could be > 3 times larger if the peak energy E p > 1-2 McV. 
Figure 9 shows the isotropic equivalent luminosity of GRB X-ray afterglows scaled to t = 10 hr after the 
burst (in the cosmological rest frame of the source), Lx(10hr), as a function of their isotropic gamma-ray 
energy release, £ 7 , iso , for GRB 050509b together with a sample of long GRBs. L x (10hr) for GRB 050509b 
is estimated by extrapolating the flux measured by the Swift XRT using the best-fit power-law decay index 
of a = 1.3, which is also consistent with the Chandra upper limit. 

A linear relation, Lx(10hr) oc -E 7j i so , seems to be broadly consistent with the data, probably suggesting 
a roughly universal efficiency for converting kinetic energy into gamma rays in the prompt emission for both 
short and long GRBs. This "universal" efficiency is also likely to be high (i.e., the remaining kinetic energy 
is comparable to, or even smaller than, that which was dissipated and radiated in the prompt emission). If 
this is the case, the well-known efficiency problem for long GRBs also persists for short GRBs. 

The X-ray luminosity at 10 hr is used as an approximate estimator for the energy in the afterglow 
shock, since (a) at 10 hr the X-ray band is typically above both v m and v c so that the flux has a very weak 
dependence on eb [to the power of (p — 2)/4] and no dependence on the external density, both of which have 
relatively large uncertainties (Freedman & Waxman 2001; Piran et al. 2001; Berger et al. 2003); and (b) at 
10 hr the Lorentz factor of the afterglow shock is sufficiently small (r w 10) so that a large fraction of the 
jet is visible (out to an angle of ~ m 0.1 rad around the line of sight) and local inhomogeneities on 
small angular scales are averaged out. Furthermore, the fact that the ratio of Lx(10hr) and E ly i so is fairly 
constant for most GRBs suggests that both can serve as reasonable measures of the isotropic equivalent 
energy content of the ejected outflow. A possible caveat to the above statement arises if the observer is 
in fact not within the aperture of the GRB jet (as is suggested to be the case in both X-ray flashes and 
X-ray rich GRBs; Granot et al. 2005). In this case i? 7 .iso can be significantly smaller than the isotropic 
equivalent kinetic energy in the afterglow shock, which is better reflected by Lx(10 hr). This is likely to 
be the reason why GRB 031203 is above the correlation shown in Figure 9 (Ramirez- Ruiz et al. 2005). An 
off-axis interpretation for GRB 050509b, on the other hand, is unlikely since its X-ray afterglow light curve 
was observed to decay from a very early epoch (at t w 10 2 s). This is also consistent with the fact that 
GRB 050509b falls close to the correlation. 

The above arguments suggest that the energy in the outflow ejected by GRB 050509b was ~ E~y [ so ~ 
10 48 - 5 erg, if it was spherical. On the other hand, if it was collimated into a narrow jet of half-opening angle 
9o, then the true energy would be smaller by a factor of / b = (1 — cos# ) ~ #o/2- Since a significant off-axis 
viewing angle is not likely, the true energy probably does not exceed E 7 ^ so . A higher redshift would increase 
2? 7l i so and with it the estimate for the energy release in this event; however, it would still remain significantly 
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less energetic than typical long GRBs (see Figure 9). 

As also argued by Lee et al. (2005), the fact that the X-ray afterglow luminosity of GRB 050509b is 
much smaller than that of long GRBs is probably because the event was sub-energetic, rather than due 
to differences on the values of the external density or the microphysical parameters. This is illustrated in 
Figure 8 by a fit to the currently available afterglow data using parameter values that are typical for long 
GRBs, except for the isotropic equivalent energy in the afterglow shock, Ek.iso: which is here taken to be 
equal to i£ 7 ; so assuming z = 0.2248. Other parameter values could also give a reasonable description of 
the rather sparse data. In Table 3 we demonstrate a few different sets of parameters that fit the afterglow 
data. Again, we refer the reader to Lee et al. (2005) for a more detailed description. Regardless of the 
redshift, it will be very difficult to detect the afterglow in the radio, since the maximal flux density (given 
the observational constraints) is unlikely to exceed ~ 15 /iJy. 

If short GRBs occur significantly outside of their host galaxies, as may be common for binary mergers 
(Tutukov & Yungclson 1994; Bloom et al. 1999; Fryer et al. 1999; Bulik et al. 1999; Belczyhski et al. 2000), 
then one might expect the external density encountered by the afterglow shock of some GRBs to be very 
low, typical of the intergalactic medium (IGM), tiigm ~ 10 _6 5 (1 + z) 3 cm~ 3 . This may help explain why 
some short bursts could have very faint afterglows. Since GRB 050509b happened to occur near the center 
of a galaxy cluster where the external density is relatively high, its X-ray afterglow was relatively brighter. 
If indeed GRB 050509b is associated with the galaxy cluster at z w 0.22, then one might expect the external 
density to be intermediate between the IGM and ISM: the GRB is ~ 76" from the center of the cluster 
as determined by the X-ray position (§ 3.4), corresponding to <~ 270 kpc in projection and well within the 
diffuse emission from the hot intracluster medium gas (which extends to a radius of ~ 1 Mpc). This suggests 
an ambient density near the position of the GRB of n w 10 _3 -10 -2 cm~ 3 , though this estimate is uncertain 
because the space position of the burst relative to the cluster center and the intracluster medium (ICM) 
density profile arc not known precisely. 



7. Discussion 

The lack of a strong afterglow signature sets GRB 050509b apart from most other GRBs. 8 As a com- 
parison, the low-redshift long-soft burst (GRB 030329, z = 0.1685; Grciner et al. 2003), if placed at the 
redshift of Gl, would have been R w 14 mag at t = 8000 s; this is approximate 10 mag brighter than the 
detection limits found herein. Even at z = 3, the optical afterglow of a GRB 030329-like burst should have 
been detected at early times (neglecting the effects of dust extinction). Our non-detections (i? ;> 24 mag) 
of variability at 1.3 hr in the what would be the restframe at z=l is more than 3.5 magnitudes deeper than 
the faintest optical tranisent found for a long GRB (021211; z=1.0; see fig. 2 of Fox et al. 2003). 

The lack of detectable optical/infrared afterglow is not surprising on grounds related to the progenitors 
and to GRB afterglow theory. First, since the luminosity of long- wavelength afterglows scales with the 
square root of the ambient density (Begelman et al. 1993; Meszaros et al. 1998), events that occur in the 
ISM or IGM should be intrinsically fainter (at optical/infrared wavelengths) than those occurring in the 
circumburst environments of collapsars (see Panaitescu et al. 2001). Second, based on < V/V max > studies, 
the isotropic-equivalent peak luminosity (L p (j)) of short bursts is similar to that of long bursts (Schmidt 



8 We note that two other low-luminosity, low-redshift GRBs (980425 and 031203) had detectable X-ray afterglow emission 
but no optical afterglow. The only transient optical signatures for these bursts were apparent supernovae. 
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2001), implying that the total energy output (-E 7j i so ~ L p (7)/r?x duration, with 77 as the conversion efficiency 
to gamma rays) is at least an order of magnitude smaller for short bursts. As argued by Panaitescu et al. 
(2001), since afterglow brightness scales with Ej^ so (l — -q), short-burst afterglows would be systematically 
faint. 

Now that there is a detected X-ray afterglow we are in a position to directly test the faintness claim, 
by inferring the gamma-ray energy release and X-ray afterglow luminosity (a proxy for the kinetic energy 
in the blast wave). From Figure 9 it is clear that this ratio for GRB 050509b is similar to that found in 
long-duration GRBs. This is a striking observational bridge to long-duration bursts and suggests a common 
physical mechanisms for prompt and delayed (afterglow) emission for both long-duration and short-duration 
GRBs, even though their progenitors are probably different. 

A tentative detection of an afterglow signal by adding up the emission of 76 short BATSE bursts was 
reported by Lazzati et al. (2001) (see also Connaughton 2002). The signal peaked at t k, 30 s after the 
burst trigger with a relatively flat vF v «5x 10 -10 erg cm -2 s . This would correspond to an X-ray flux 
in the 0.2-10 keV range, of F x « 2 x 10 -9 erg cm" 2 s" 1 . The X-ray flux of the afterglow of GRB 050509b 
is best constrained around t s=s 200 s, and is found to be Fx ~ 6.5 x 10 -13 erg cm -2 s _1 . Extrapolating 
this flux to t w 30 s with a power-law index in the range inferred from the data, 1.0 < a < 1.7, gives a flux 
that is lower than the one found by Lazzati et al. by a factor of <~ (1 — 5) x 10 2 . This might suggest that 
cither the possible detection by Lazzati et al. (2001) was not statistically significant, or the X-ray afterglow 
of GRB 050509b is underluminous compared to the average value for short GRBs by at least two orders of 
magnitude. 

With essentially no indication of recent star formation in Gl, massive progenitor stars leading to col- 
lapsars cannot be present in Gl. SI and S2, the brightest and third-brightest sources within the XRT error 
circle, have no indication of recent star formation if their rcdshifts arc < 1.3 (SFR < 0.05M Q year -1 for 
z < 0.3 and SFR < 1M© year -1 for z < 1.2). The fainter (and blue) objects discussed in § 5.1 are likely to 
be background galaxies. If the origin of GRB 050509b is from a collapsar, it is likely that its redshift exceeds 
1.3. 

If GRB 050509b is a background object at z > 2, some progenitor scenarios are difficult to reconcile. 
With an observed duration of <~ 30 msec, the rest-frame duration would be only about 10 msec. This is 
implausibly short for an NS-NS merger, and marginally possible for a BH-NS merger if the coalescence is 
through unstable mass transfer (Lee & Kluzniak 1999; Rosswog 2005; Miller 2005). It is hard to simulta- 
neously accommodate the short intrinsic timescale and the higher energy budget of the burst within any 
compact merger model, if it is at high redshift. 

If short GRBs trace star formation with a time delay through double compact mergers with coalescence 
time scales of 10 7 -10 10 yr (as opposed to prompt tracers of star formation as with the collapsar scenario 
for long GRBs; Bagot et al. 1998 and Bloom et al. 1999), then we expect some fraction (10-30%) of short 
GRBs to be seen in association with early-type galaxies in general and clusters specifically (see Nutzman 
et al. 2004 for rate density in the local universe). This is somewhat model dependent, since the distribution 
of compact merger timescales is poorly constrained by data, but broadly consistent with both observed and 
model distributions. 

A core-collapse supernova (SN) produces no electromagnetic radiation until its envelope is completely 
consumed by the explosion (although see Khokhlov et al. 1999). This phase ends, however, with a brilliant 
flash of X-ray or extreme ultraviolet photons as the shock reaches the stellar surface. The "breakout" flash 
is delayed in time, and vastly reduced in energy, relative to the neutrino transient produced by core collapse. 
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However, it conveys useful information about the explosion. Shock breakout flashes were predicted by Colgate 
(1968) as a source for (the then undetected) gamma-ray bursts. The explosion of SN 1987A stimulated a 
reanalysis of supernova breakout flashes by Ensman & Burrows (1992) and, more recently, by Blinnikov 
et al. (1998) and Blinnikov et al. (2000). These studies represent an increase in sophistication toward the 
full numerical treatment of this complicated, radiation-hydrodynamic problem. In principle, the XRT data 
could constrain the existence of a shock breakout produced by both a red supergiant explosion like SN 1993J 
(Van Dyk et al. 2002) and a blue supergiant explosion analog to SN 1987A, but the X-ray luminosity is 
sensitive to the uncertain distribution of extragalactic gas column and the specific XRT observing epochs. 

Using our ESI optical imaging, we can also limit the presence of brightening due to a supernova or 
supernova-like emission at 8.17 d after the GRB to R c ss 25.0 mag. A normal, unextinguished Type la 
(thermonuclear) supernova at z — 0.22 would have R 22 mag, around 6.7 d after explosion (t = 8.17 d 
in the observer's frame). A very subluminous SN la like SN 1991bg (Filippenko et al. 1992) would have 
R s=a 24 mag, still somewhat brighter than our limit. Extinction would obviously make the SN fainter, but 
the Milky Way contribution is small (Ay ~ 0.06 mag; Schlegel et al. 1998), and the outskirts of an elliptical 
galaxy in a cluster should have essentially no dust. While some core-collapse supernovae could be as faint 
as (or fainter than) our limit, the presence of such a supernova in the outskirts of an elliptical galaxy would 
be truly extraordinary (see van den Bergh et al. 2005). Others have also reported no evidence for a SN at 
later-times (Hjorth et al. 2005; Bersier et al. 2005). 

The location of this (and future) short burst provides a useful discriminant for distinguishing between 
different progenitor models of short bursts. Simplistically, we would expect evaporating black holes to occur 
near the center of deep potential wells (as discussed in the context of Galactic BHs; Cline et al. 1999); thus, 
the offset from Gl seems to disfavor this hypothesis. A giant flare from a magnetar would need to have a 
isotropic luminosity (£ 7 .i so ) larger by a factor of ~ 10 3 and an -E 7j i so larger by a factor of ~ 10 2 compared the 
the initial spike of the 27 December 2004 giant flare from SGR 1806-20 (the difference in the factor between 
the two quantities arises since GRB 050509b lasted only ~ 30 ms, which is <~ 10 times shorter than the 
initial spike of the giant flare from SGR 1806-20). Bursts from magnetars might be expected from later-type 
galaxies than Gl where neutron stars would be formed copiously: magnetic field decay would cut the active 
lifetime for megaflare activity after <~ 10 4 yr. 

8. Conclusions 

We have monitored the location of GRB 050509b at optical and infrared wavelengths from 8 minutes 
to 8 days after the trigger and found no indication of variability at the location of the fading X-ray source, 
the first solid X-ray detection of an afterglow of a short-hard burst. Near the location of this source we and 
others have found an apparent group of faint blue galaxies at redshifts > 1.3. While it is indeed plausible 
that this short burst arose from a progenitor connected with those galaxies, we found — based on a positional 
argument — plausible evidence that the progenitor is likely associated with Gl, a bright elliptical galaxy 
at z = 0.2248. We have argued that the observations find natural explanation with a compact merger 
system progenitor. If so, then short-hard GRBs provide a bridge from electromagnetic to gravitational wave 
astronomy: indeed has GRB 050509b occurred a factor of ^3 closer in luminosity distance is might have 
produced a detectable chirp signal with the next generation LIGO-II gravitational wave facility 9 . 



9 http://www.ligo.caltech.edU/docs/G/G9901 1 l-00.pdf 
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Brightcning emission from most types of supernovae would have been seen in our imaging, so the lack of 
such emission appears inconsistent with the notion that short bursts are due to collapsars or variants thereof. 
Our afterglow modeling is also consistent with, though does not require, a circumburst medium having lower 
density than that inferred in long-duration GRBs; if true, this would suggest that the progenitor produces 
a GRB in an environment that is baryon poor compared to that expected for collapsars. Moreover, we have 
seen no evidence for ongoing star formation in the putative host, so there are likely no remaining massive 
stars. Given the short active life of a neutron star having a high magnetic field, this also disfavors the 
magnetar hypothesis. 

The non-detection of brightening emission may place limits on the presence of a thermal "mini-supernova" 
from non-relativistic ejecta of a compact merger system (Li & Paczyhski 1998; Rosswog & Ramirez-Ruiz 
2002). In this scenario, the small dense mass (m e j) ejected during coalescence expands as it is heated by ra- 
dioactivity of the decompressed ejecta. Using the scalings of Li & Paczyhski (1998) and crudely assuming that 
10% of the bolometric light at peak is radiated in the R band, the i?-band brightness should peak at observer 
time t w l^TOoj/O.OlMo) 1 / 2 d after the burst, with absolute magnitude M R w -18.5-1.25 log(m cj /0.01M o ) 
mag. Assuming that the GRB did indeed originate from the redshift z = 0.2248, upon inspection of Figure 
1, with non-detections at Mr w —16 mag at t w 1 d, we can very roughly exclude m c j >few xl0~ 3 M Q . 
Though the Li & Paczyhski (1998) model was intended as a simplistic sketch of the phenomenon, this limit 
on 77i j is somewhat surprising given the amount of escaping non-relativistic material expected in compact 
mergers (Rosswog & Ramirez- Ruiz 2002). Indeed, we consider this lack of a "mini-supernova" as weak evi- 
dence against a z — 0.22 origin from a compact merger system. Still, these limits are subject to considerable 
uncertainty in a number of uncertain parameters of ejecta. For instance, if the velocity of the ejecta were to 
be ~0.01c instead of 0.3c (as assumed by Li & Paczyhski 1998) then the peak of the thermal emission would 
occur after about 1 month, and would not have been detected with the current limits. 

We conclude by emphasizing that in the NS-NS or BH-NS progenitor hypothesis for short-hard bursts, 
the hosts galaxies may be a range of Hubble types (e.g., Livio et al. 1998). Compact merger systems coalesce 
in appreciable rates from Myr to Gyr after a starburst (e.g., Fryer et al. 1999; Bloom et al. 1999). Obviously, 
the longer the time since the starburst, the larger the distance a binary system will travel before coalescence. 
A clear prediction from this model is that as more short bursts are localized, those associated with later-type 
galaxies of a given mass should be preferentially closer to the star-formation centers of the host; that is, we 
expect a more concentrated distribution around a spiral galaxy with the same mass as an early-type. On 
the other hand, dwarf star-forming hosts have shallow enough potentials that merger systems from these 
galaxies could coalesce at appreciable distances (;> lOOkpc) even shortly after starburst. As Swift localizes 
more short-hard bursts, we expect that the offset distribution around galaxies will further elucidate the 
progenitor question. 
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Fig. 1. — Upper limits from WIYN, PAIRITEL, and Keck images (filled symbols) along with those reported in 
the literature (open symbols). Magnitudes are corrected for extinction using the maps of Schlegel et al. (1998) 
and converted to AB magnitudes, assuming z = 0.2248 (see text), using the relations of Frei & Gunn (1994) and 
Blanton et al. (2003a). Times are reported in seconds relative to the Swift trigger on May 9.166889. The WIYN 
and PAIRITEL upper limits are 5a, but many of the quoted limits in the literature are not accompanied by a stated 
significance level. For smaller telescopes with large pixels, the light from the nearby galaxy is likely a significant 
contaminant, resulting in upper limits that may be overestimated in the literature. 
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Fig. 2. — Keck/LRIS G-band (top) and _R-band (bottom) images, zoomed to show the XRT error circle. The larger, 
blue circle is the revised XRT position from Rol et al. (2005); the green circle to the west and north of that is our 
2a confidence region (see text) of the XRT position. The images have been median subtracted and smoothed to 
accentuate the detection of faint sources under the glare of the bright galaxy Gl. The 11 sources consistent with 
the Rol et al. (2005) X-ray afterglow localization are labeled in both images, with the astrometric positions given in 
Table 2. North is up and east is to the left. Gl is the large galaxy to the west and south of the XRT. Bad pixel 
locations are denoted with "BP." 
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Fig. 3. — Scatter plot of the position offsets between the XRT sources (other than the afterglow candidate) 
and their nearest sources in the wide-field Bok image B-band image. Two of the XRT sources each have two 
Bok sources nearby. For these two XRT sources, the dashed line represents the offset to the closest source, 
and the dotted line represents the offset to the next closest source. There is strong clustering around an 
offset of 4.5" in RA indicating a global shift in the absolute astrometric zero-point of the XRT frame. The 
number of X-ray/optical cross-matches is reasonable given the sensitivity of the two frames (see text). 
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Fig. 4. — Compilation of observed X-ray light curves for GRBs and core-collapse supernovae (adapted from Kouve- 
liotou et al. 2004). The X-ray luminosities for the short GRB 050509b are calculated assuming z = 0.2248. For GRB 
050509b we used data points (square symbols) from both the XRT and the Chandra upper limit reported by Patel 
et al. (2005). The solid line corresponds to the best-fit temporal decay index of —1.3. 
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Fig. 5. — A 14' x 14' adaptively smoothed image of the XRT data. GRB 050509b is embedded in the diffuse 
X-ray emission associated with the galaxy cluster, and is about 270 kpc in projection from the cluster center. 
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Fig. 6. — Keck/DEIMOS spectrum of the galaxy Gl (along with its variance spectrum) located ~ 10" west 
of the center of the XRT error circle for GRB 050509b. The data were obtained using the 600 line mm -1 
grating centered at 7200 A and the galaxy was observed through a 1.1" slit (FWHM rs 5 A). The strong 
absorption-line features indicate z — 0.2248, and a comparison of the spectrum against a template spectrum 
of HD 72324 provides an estimate of the velocity dispersion: a = 260 ± 40 km s _1 . 
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Fig. 7. — False-color image of the field of XRT constructed with the G (blue) and R (red) Keck/LRIS images; 
green is interpolated between the observed bands. Aside from source S6 (which appears red) and possibly 
the J2/J4 complex, all of the XRT-consistent sources appear to be faint and blue, consistent with a small 
group of star-forming galaxies at a redshift larger than the cluster. As seen in this image, a number of such 
groups appear throughout the field (there are 2 faint blue galaxies to the north of Gl, also embedded in the 
light of Gl). 
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Fig. 8. — The spectrum at t = 80 s after the trigger (upper panel) and the light curves (lower panel) of the X-ray 
(1 keV; blue solid line), optical (R band; red dashed-dotted line), and radio (4.9 GHz; green dashed line) emission for 
a spherical afterglow shock propagating into a uniform ambient medium, using the model of Granot & Sari (2002). 
We also show the flux normalization in the X-rays from the Swift XRT detection, as well as upper limits in the 
optical and in the radio. Here, we adopt the redshift of the tentative host galaxy (z = 0.2248) with typical interstellar 
medium density n = 1 cm" 3 , and assume that the isotropic equivalent kinetic energy in the afterglow shock is equal to 
B T ,iso (i.e., 2.7 x 10 48 erg). The microphysical parameters are taken to be typical of those inferred from the modeling 
of afterglows of long GRBs: p = 2.2, es = 0.15, es = 0.046. Table 3 gives other models consistent with the data. 
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Fig. 9. — Isotropic equivalent luminosity of GRB X-ray afterglows scaled to t = 10 hr (source frame) after the burst 
as a function of their isotropic gamma-ray energy release (adapted from Kouveliotou et al. 2004) . If GRB 050509b is 
located at z = 0.2248, the isotropic equivalent luminosity of the X-ray transient at t = 10 hr assuming t~ 13 and the 
isotropic gamma-ray energy would be ~ 9 x 10 40 erg s _1 and ~ 2.7 x 10 48 erg, respectively (black symbol). 
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Tabic 1. Upper Limits on Optical-IR Afterglow of GRB 050509b 
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a Time since 04:00:19 UT, the time of the Swift/BAT trigger. 

b imits were converted to AB magnitudes following (Frei & Gunn 1994). Limits 
are presumed to be 5a detection limits for a point source (not all reports in the 
literature stated the significance of the upper limits). No extinction correction 
has been applied to these magnitudes. 
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Table 2. Properties of Faint Sources in the XRT Error Circle 
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a J2000; Positions are based on image centroids and have been tied to the SDSS astrometric frame 
(see text). The absolute positional uncertainty, inherited in part from the SDSS uncertainties, is 
0.2". The relative uncertainty between object positions is expected to be ~0.1". 

b Photometry has been tied to the SDSS calibrations of the field. We used a color term in the 
conversion of the Keck R to the Sloan r' filter and constant magnitude offset (following Windhorst 
et al. 1991). We assumed no color term in converting G to g' . Errors include the 0.1 mag zero-point 
uncertainty in the tie to the SDSS photometry. We do not report photometric errors for those sources 
near the detection limit, but such errors are likely to be of order 50%. No extinction correction has 
been applied to these magnitudes. 
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Table 3. Representative Fits to the Afterglow of GRB 050509b 



z 


E ktiso /erg 


n/cm 3 
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Note. — An example of four different sets of parameters that 
fit the rather sparse afterglow data for GRB 050509b. In the first 
two cases we fix the redshift at that of the putative host galaxy 
(Gl; z = 0.2248), and assume an external density that is either 
typical of the interstellar medium (ISM), 1 cm~ 3 , or typical of the 
intergalactic medium (IGM), 10~ 6 cm~ 3 . In the last two cases we 
explore the option of a relatively high redshift, z — 3 with the same 
two very different values for the external density. In all cases, the 
values of the micro-physical parameters (p, e e and e B ) were chosen 
to be typical of those inferred from afterglow fits for long GRBs. 



